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We pick up a method originally developed by Cheng and Tsai for vacuum perturbation theory 
which allows to test the consistency of different sets of Feynman rules on a purely diagrammatic 
level, making explicit loop calculations superfluous. We generalize it to perturbative calculations in 
thermal field theory and we show that it can be adapted to check the gauge invariance of physical 
quantities calculated in improved perturbation schemes. Specifically, we extend this diagrammatic 
technique to a simple resummation scheme in imaginary time perturbation theory. As an application, 
we check up to 0{g'^) in general covariant gauge the gauge invariance of the result for the QCD 
partition function which was recently obtained in Feynman gauge. 
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I. INTRODUCTION 

In finite temperature quantum field theory for massless degrees of freedom, the loop expansion does not coincide 
with an expansion in orders of the coupling constant. Naive perturbative calculations lead to gauge dependent and 
infrared divergent results for physical quantities [1]. This necessitates the reorganization of the perturbative expansion 
in consistent resummation schemes. One of the most important consistency checks for such schemes is the test of their 
gauge invariance. In practice, this is done mostly by explicit calculations [2], though Ward identities can be used as 
well [3]. In the vacuum sector, yet another technical tool is at hand. This is a purely diagrammatic method developed 
by Cheng and Tsai [4] which allows to establish the gauge invariance of a set of Feynman diagrams without carrying 
out explicit loop calculations. Recently a similar diagrammatic method was independently developed by Feng and 
Lam [5]. The main motivation of the present work is to undertake a first step in extending the Cheng/Tsai method 
to resummed perturbative calculations. To this end, we review and slightly extend the approach of Cheng and Tsai 
presenting a complete set of diagrammatic rules. As an example we show how to check in this approach the gauge 
invariance of the two loop contribution to the vacuum partition function to 0{g'^) in general covariant a-gauge. Then, 
we extend these rules to a simple resummation scheme. 

As a first application, we focus on the gauge invariance up to 0{g^) of the QCD partition function. The latter 
has recently been calculated up to 0{g^) [6-8] using a gluon propagator with a resummed static mode. Indeed, there 
are several reasons for choosing this example. First, the calculations to order 0{g^) and 0{g^) exist in Feynman 
gauge only, and the authors did not check the gauge invariance of their results. Since the chosen resummation scheme 
amounts to the "free" Lagrangian of a massive Yang-Mills field, which is not gauge invariant, invariance of the fi- 
nal result is not automatic. This is different from naive perturbation theory where gauge invariance of observable 
quantities is automatic if one uses a complete and consistent set of Feynman rules. (For certain gauges this may 
be a non-trivial problem in itself.) In fact, the method of Cheng and Tsai was invented to study the consistency of 
different sets of Feynman rules in this respect. Second, in contrast to dynamic quantities of finite temperature QCD 
which in general require resummed propagators and vertices [2] , the calculation of the partition function requires only 
resummed propagators. So, our work extends the method of Cheng and Tsai to the simplest case of a resummation 
scheme. 

In Sec. II we outline the method of Cheng and Tsai [4] and give a short example. In Sec. Ill we discuss the resum- 
mation scheme and some general aspects of the free energy (resp. partition function). As an application we extend 
the example discussed in Sec. II to finite temperature and prove the gauge invariance of the QCD free energy up to 
0{g^). In Sec. IV we discuss a technical issue which arose in the Sees. II and III, the problem of shifting momentum 
variables. In Sec. V we then prove the gauge invariance of the free energy in finite temperature QCD up to 0{g^). 
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II. THE DIAGRAMMATIC METHOD OF CHENG AND TSAI 



In the Feynman rules of QCD, the chosen gauge manifests itself in the Z;:^ -dependent part of the free gluon propagator 
D'^Jl and in the way D'^Jl couples to ghosts. We write the gluon propagator in the general form 



k'^ + is 



= Df(k) + Al''{-k)K-At\k)k,. 



(2.1) 



+ 



Diagrammatically we denote the factor k^ by an arrow-head. The Feynman-propagator D"^^^ is given by 



k'^ + is 



and Af,^ is defined as 

fj, 



(2.2) 



(2.3) 



The gauge dependent pieces are contained in the function ai^{k), e.g. a^(/i:) = ^{1 — a)^ for covariant a-gauges. 
Checking the gauge invariance of a set of Feynman graphs amounts then to the assertion that the result is unaffected 
by changing a^(k), i.e. by choosing a different gauge. In the following we review the method of Cheng and Tsai [4] to 
do this check diagrammatically. 

The free gluon propagator in (2.1) has been separated diagrammatically into its Feynman gauge part Dp and 
explicitly gauge dependent parts. In this way the method of Cheng and Tsai allows to disentangle arbitrary diagrams 
into their Feynman gauge parts and remainder terms. A set of diagrammatic rules is used to test whether the sum 
of these remainder terms vanishes. If this is the case, the calculation is independent of the particular choice of (/;:). 
The original motivation of Cheng and Tsai was to check by this technique the consistency of the Feynman rules for 
different gauges [4,9,10]. We shall employ the same technique to show the gauge invariance for quantities calculated 
in resummation schemes. 

A. The Basic Idea 

We start with —iS^^ ^k^^ie' ^ gauge-dependent part of the propagator. We split this term into the two factors 
A^^{k) and kj,. If this is a contribution to an internal line of a Feynman diagram, then kj, is Lorentz-contracted with 
a vertex. For the case of a 3-gluon-vertex, this contraction results in a very simple structure. 



(2.4) 



Each one of the two terms in square brackets has the structure of a transverse projector e.g. Vi^p{q) = Qiip—Qii^p^ 
acting on one of the 2 outgoing lines of the vertex. Note that in what follows we refer to P as a "projector" though 
it does not satisfy the normalization, i^V . Diagrammatically, we express (2.4) in the following way: 




+ 





(2.5) 
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The blob at the vertex on the l.h.s. indicates that this vertex carries a complete Lorentz structure which has 
disappeared on the r.h.s. upon contraction. Note that the diagrams on the r.h.s. have opposite sign. 
V acts on the next propagator in the following way: 

gr'"'P'^{p)Df,{p) = - igr'^g^, - Gf%p)p,. (2.6) 

We represent this equation graphically as 




The first term on the r.h.s. of (2.6) is momentum independent; sandwiched between two vertices attached to either 
end of the line it leads to a contraction of the vertices (cf. (2.17)). It is graphically represented by the so called 
"contraction arrow". The G in the second term on the r.h.s. is nothing but the product of a ghost-gluon vertex and 
a ghost propagator 

Gf\k) = -igr"^ [(a-k- 1) - fc^a^] (2.8) 
and is therefore called the ghost contribution. For any covariant gauge G{k) reduces to 

GZ(.k)=igr''^^. (2.9) 
Putting (2.4) and (2.6), respectively (2.5) and (2.7), together, we have 




(2.10) 



There are three more diagrammatic rules needed for the decomposition of Feynman diagrams into explicitly gauge 
dependent and gauge independent parts. The first specifies how V acts on A: 

gr'>'Vj;{q)A';f{q) = Gf\q) - Gfj{q). (2.11) 

Note that for all covariant gauges Gf = G and thus the r.h.s. is zero. Therefore we get for covariant gauges using 
(2.4), (2.6) and (2.11) 




(2.12) 



We will work in general covariant gauge in the following applications, so rule (2.12) will be often used. 
The second rule concerns the case that an arrow-head k/^ acts on a ghost-contribution. It reads 

rGf'^(^) = -^''Gf<=(g) -p''Gf<=(g) 

= -igr''-p''Gt''{q). (2.13) 

Graphically we have 
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(2.14) 



where we denote igf'-^'"^ by 1. This relation is a direct consequence of energy-momentum conservation at the vertex. 
The last rule concerns the case that is connected to a fermion-gluon vertex F. We have 



= igin 



- rrif -\- IS 
1 

— rrif + is ']/) — rrif -\- is 



(jl)-mf ^ is) 
1 



which win be represented graphicahy as 






(2.15) 



(2.16) 



In the last step of eq. (2.15), we have used k + q = p. 



B. The Basic Identities 



In this subsection we give rules which allow to relate gauge dependent parts of diagrams with different topology 
with each other. This will be used to show that the gauge dependent parts of different diagrams cancel each other. 
The first rule relates a diagram with a 4-gluon- vertex Q to those with contraction arrows. It is obtained by applying 
(2.10) to diagrams with a second three-gluon vertex F: 



/ia 



z/b k 





Sd fidi Sd /ia 



e ~-k' z/b k ~4>e -k' 7C j/b k -k' 7C j/b k 




6d 
= 0. 

k' 7C 



i99''^r'mJ-p - k + k',p + k, -k') 



-igg-><t>r-rZ%i-p -k + k', k,p - k') 

-igg''f'f^-'ri%{k' -k,k,-k') 

The second identity addresses the case that the contraction arrow points at a four-gluon vertex: 

/ia g Sdfidi Sdfidi Sd /ia Sd 



= 0. 






7Cz/b 



7Cz/b 



7C z/b 



7C 

fbhg ^u(f) ^adch q 



(2.17) 



igr'^ag^'fQhdcb^ - ig^'^g'^QX, - igf^^g'^^Qltl + i9r'9''^QX'i = ^- (2-18) 

The third identity (which we have not found in the literature, but which will play an important role in Sec. V) 
concerns the case that the contraction arrow points at a vertex which carries a G. This identity follows from the 
Bianchi identity for the color structure coefficients f^^^: 
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CfcGi d c (j)e d c^^G d 



igr'Gl^'ik) - i9f''G'^''"'{k)g,^ - i9r'"G'^'''^\k)g,^ =0. (2.19) 
The fourth identity treats the case that the contraction arrow points at a fermion-gluon vertex: 



a i a i a i 




-i9r''g'"^Fl,j{p + k)= Tl% - Tl%. (2.20) 

Note that the gluon hnes at the four-point fermion-gluon vertex created in (2.16) cannot be interchanged due to 
the color-structure (cf. Appendix C). 

C. A First Example 

To illustrate the method of Cheng and Tsai, we discuss briefly the gauge invariance of the 0{g^)^ contributions to 
the QCD partition function. For this purpose we will from now on concentrate on general covariant a -gauge. That 
means that we will use the special structure Gf = which the ghost contribution satisfles (cf. 2.11). This leads to 
the use of eq. (2.12) in the following. The partition function consists of four diagrams ^: 




(2.21) 



The dotted lines are the usual notation for ghost propagators. For our purpose it is however more useful to express 
diagrams like (2.21.C) in the following way: 




(2.22) 



We follow a three step procedure: 
We select, one after the other, a diagram in (2.21) and separate the a^-dependent part of one of the internal gluon 
lines. E.g. for diagram (2.21.B) this looks like 



^We indexed each diagram with a letter in order to refer to it later in the text (e.g. diagram (2.21.C) is the third diagram in 
(2.21)). Moreover we give the factor (-1) for the fermion and ghost loops in the following explicitly in the symmetry factor. 
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e 



D 

(A) 




+ 





(2.23) 



(B) 



(C) 



(D) 



Because A is proportional to an odd power of the momentum k diagrams (2.23.C) and (2.23.D) can be summed. 
So we get: 



e 



D 

(A) 




+ 2 




(B) 



(C) 



(2.24) 



Next we use the rules of Sec. II A to decompose diagram (2.24.C): 




(A) 



(2.10) 



+ 




+ 






+2 





(2.10) 

+ 2 





e 



^'^'^ + 2 




(2.25) 



As a second step, we invoke the identities of Sec. II B as follows. Connecting the external lines of eq. (2.17), we get 
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= 0. (2.26) 



(A) (B) (C) (D) 

First we observe that diagram (2.26.C) is zero due to the color structure. Obviously diagram (2. 26. A) and (2.26.B) 
look very similar, and after substituting —k to k in diagram (2.26.B) we arrive at 




(2.27) 



Clearly momentum shifts do not affect the value of a complete Feynman diagram. Here, however the integrand 
of a complete Feynman diagram is rewritten algebraically as a sum and the sequence of the summation and the 
integration over the internal momentum loop variables is interchanged. Performing the same change of integration 
variables in all terms of the sum does not effect the result. However, changing integration variables only in some 
terms, as e.g. in eq. (2.27), requires a careful study of the corresponding terms in the sum. Here, this amounts to the 
demand that a^(/i:) is not too singular to allow for this substitutions and shifts (c.f. [4]). In section IV, we shall show 
that this subtlety does not obstruct the arguments made in what follows. Altogether we get from (2.26) using (2.27) 
and attaching a A 




(2.28) 



(A) (B) 

Similarly, eq. (2.19) leads to the identity 




(2.29) 



(A) (B) 
This means that (2.25.C) cancels (2.25.D). Using (2.28), eq. (2.25) is thus simplified to 




(2.30) 



(A) (B) (C) 
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Obviously this relates the three different diagrams (2. 21. A), (2.21.B) and (2.21.C). 

We are now still left with the diagram (2.21.D). Using (2.15) we can make the following decomposition (2.16): 




(2.31) 



(A) (B) (C) 

Closing identity (2.20) on the right and left side to loops and attaching a A one can see that the r.h.s. of eq. (2.31) 
is zero and in this way the a^-dependent part of (2.21.D) vanishes. 

In the third and final step, one just has to count symmetry factors in order to show that all a^-dependent terms in 
(2.21) cancel each other. To do so we count the different ways to obtain the diagrams in (2.30) from (2.21). There are 
6 possibilities to assign A/^kjy to diagram (2.21.B): there are 3 lines to choose from and 2 ends of the line to which the 
factor A can be attached. (The latter factor 2 is explicitly given in (2.24).) Similarly, there are 4 such possibilities 
for diagram (2. 21. A) and 2 for (2.21.C). This counting shows immediately that all terms linear in A/^kjy cancel each 
other. Similarly one can show the same for the terms which are quadratic and cubic in A/^kj^. That is all we have to 
show: (2.21) leads to the same result whether we calculate it in an arbitrary a-gauge or in the Feynman gauge. 



D. The Symmetrization Argument 

Arguing separately for each order of A/^kjy is very hard to do for more complicated diagram. In this subsection we 
outline a technique for solving this problem which is based on the following observation: In principle gauge invariance 
of a sum of diagrams requires only that this sum does not change if is simultaneously changed in all lines of the 
diagrams. However, we will show in the following that the same sum of diagrams remains unchanged even if is 
changed in only one line in each diagram. To see this we start again with the sum given in (2.21) and symmetrize the 
diagrams with respect to the gluon propagators: 




(E) (F) (G) 
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- 2 



V 




(2.32) 



In the spirit of Sec. II C, we choose now one hne, say Di, and spht it into its gauge dependent and independent 
parts. This means e.g. for (2. 32. A): 



Clearly all diagrams in (2.32) containing Di must be treated in the same way. According to (2.31) the fermion 
diagram cancels immediately. By complete analogy with our treatment in eq. (2.25), one observes that the gauge- 
dependent part arising from the propagator Di in (2.33) cancels the corresponding gauge-dependent contributions of 
Di in (2.32.B/C/E). Consequently, we can replace in (2.32) the propagator Di by its Feynman part without 
changing D2 and D3. In the procedure the explicit form of D2 and D3 is never used. As we have changed Di to 
Di^F without changing the value of the whole sum of diagrams we can do the same for D2 and finally also for D3. 
Consequently, D2 and D3 can also be replaced by their Feynman parts D2,f and D^^f without changing the value 
of (2.32). The set of diagrams has the same value in Feynman gauge and in all a-gauges and in this sense it is 
gauge-invariant . 

From the fact that Di can be replaced by without using the explicit form of D2 and D3, we draw two 

conclusions: 

1. It is sufficient to check the invariance of a set of Feynman diagrams under the replacement of only one propagator 
in each diagram by its Feynman part. The rest can be handled by symmetrization. 

2. Since no assumption is made about the form of the D^, each could enter the symmetrized expression with a 
different A/^kjy without affecting the result of the calculation. 

In the next section we will generalize the diagrammatic rules presented here to the case of finite temperature. 



III. THE CHENG-TSAI METHOD FOR IMAGINARY TIME PERTURBATION THEORY 

So far, we have reviewed and completed the basic Cheng-Tsai rules for vacuum QCD. Now we explain in a first 
step how these rules get modified in imaginary time perturbation theory at finite temperature (i.e. with a compact 
imaginary time coordinate). Subsequently, the rules are extended to a resummed perturbation scheme in which a 
dynamically generated Debye screening mass is introduced for the static temporal gauge propagator. This paves the 
way for checking the gauge invariance of the 0{g^) result for Zqcd- 



The Feynman rules for imaginary time perturbation theory at finite temperature differ from the vacuum ones 
essentially only by changing the energy component of the four momentum k to discrete imaginary Matsubara 
frequencies, k^ = i2n7rf3~^ for bosons, ko = i(2n + l)7r/3~^ for fermions with P being the inverse temperature. This 
leads to the replacement of loop integrals in momentum space by 




(2.33) 



A. Diagrammatic Cheng-Tsai Rules at Finite Temperature 
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f d*k f (Pk 

The energy-momentum conservation at each vertex is modified accordingly: 

i{2i:Y5{k) /3{27r)^SnoS{k). (3.2) 

Up to factors i coming from the imaginary ko components, the propagators and vertices look formally the same as 
in the vacuum case, now being functions of k and ko (cf. Appendix B). They can be split up diagrammatically into 
their Feynman gauge and explicit a^-dependent parts in exactly the same way as before, e.g. 

= DF,^u + ^^i{-k)k^-K{k)k^\ with k"^ = f + kl, (3.3) 



where now the Feynman propagator is defined by 



I^f'm.=^''V^ (3-4) 



while Af,^ takes the form 



^f{k) = 8^'^. (3.5) 

With these changes, the diagrammatic rules (2.10), (2.12), (2.14) and (2.16) look the same, the ghost contribution 
reading now (cf. (2.8)) 

Gf\k) = -ir'' [{a-k- 1) k^ - k'a^] ^. (3.6) 

Also, the basic identities (2.17-2.20) carry over without further modification. Note that the application of the 
Cheng-Tsai method in imaginary time perturbation theory proceeds as in the vacuum case. It is not obstructed by 
the replacement (3.1). All rules are purely algebraic, establishing the cancellation of a^-dependent terms on the level 
of the integrand before any loop integrations are performed (for a detailed discussion, cf. Sec. IV). Consequently, they 
do not depend on changing the ko-\oop integral to an infinite sum. 

B. A Simple Resummation Scheme 

As already mentioned, naive calculations in finite temperature perturbation theory with massless degrees of freedom 
can lead to gauge dependent and infrared divergent results for physical quantities since the loop expansion does not 
coincide with an expansion in orders of the coupling constant [1]. Depending on the physical quantity one wants to 
calculate, this requires the use of more or less refined resummation schemes [2,11,12]. A prominent example is the 
QCD partition function Zqqb for which naive perturbation theory leads to an incomplete two-loop result, missing the 
0(^^)-contribution completely, and becomes infrared divergent at the three-loop level. At least up to order 0{g^), a 
consistent resummation scheme curing all infrared divergences is obtained by reorganizing perturbation theory with 
a massive static Ao-pi"opagator [7,11]. 

C=(^C + Im^Ag^^^/fco.o) - im^AgAg^o.o. (3.7) 
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Here, m denotes the dynamical generated Debye screening mass which takes into account the screening of electric 
fields in a polarizable plasma at finite temperature and density. At lowest order in the coupling constant the Debye 
mass is obtained as the infrared limit of the corresponding component of the one-loop self energy 11^^, 



m = lim Iloo{ko = 0,k). 



(3.8) 



-n 



(3.9) 



In the improved imaginary time perturbation theory, resummed according to (3.7) and (3.9), the gauge boson 
propagator reads 



ah 



m? + k'^ 



1 



F (m2 + F) 



+A^{-k)k^ - A^{k)k^. 



(3.10) 



D' 



+ 



(3.11) 



The counterterm in (3.7) is treated diagrammatically as a new two-point vertex 



5"" Si,oSi_,oSkoom 



(3.12) 



Note that the counter term in (3.7) breaks the gauge invariance explicitly. Hence, it is a priori unclear whether 
this reorganization of perturbation theory preserves the gauge invariance at a fixed order in the coupling constant. 
Due to the mass term in (3.10), the Cheng-Tsai rule for the action of the "projection" V on the propagator DJ^J, 
receives now an extra contribution 



(3.13) 



V 



+ 



G 
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Thus (2.10) is modified to 




(3.14) 



Since eq. (3.7) does not involve resummed vertices, all Cheng-Tsai rules relating different vertices to each other 
remain unchanged. Consequently, we have the diagrammatic reformulation of the energy-momentum conservation 
(2.14), the fermion- vertex rule (2.16) and the basic identities (2.17-2.20) at our disposal. Additionally, there is a rule 
linking the counter term (3.12) to the mass-dependent term in (3.13): 

2 

771 ~r 1^ 



C. The QCD Partition Function 



With the improved imaginary time perturbation scheme described in the last subsection, the QCD partition function 
has been calculated to 0{g^) [6-8]. Since it is our main aim to check the gauge invariance of these results, let us 
shortly recall them. 
Zqcd is given by the trace 

Zqcb = Tre-^^ (3.16) 
where the Hamiltonian is defined in terms of (3.7). The corresponding free energy is 

F = -|lnZQCD. (3.17) 

Its lowest order contribution is the Stefan-Boltzmann part, which is diagrammatically represented through one loop 
diagrams and has been calculated in many different gauges (c.f. [13]). The higher order corrections to this free part 
can be represented diagrammatically up to 0{g^) as follows: 
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® 



(K) 




(3.18) 



Here, the Feynman rules of section IIIB are understood, and the blob is defined in terms of the self-energy of (3.9) as 



n 



+ 



(3.19) 



The complete 0{g^) result of these contributions has been obtained by Kapusta [14] nearly two decades ago. It is 
given via the resummed two-loop calculation including the diagrams (3.18.A-E), and was obtained in [14] in Feynman 
gauge. Only recently, Arnold and Zhai have calculated the 0{g^) contributions [6,7], being able to reduce all occurring 
loop integrals to standard ones and to the so-called scalar basketball diagram. At this order all mass terms can be 
neglected in the 3-loop diagrams. With essentially the same methods but including the mass terms in the three-loop 
diagrams, the calculation was pushed to 0{g^) by Zhai and Kastening [8]. Beyond 0{g^), simple resummation based 
on (3.7) is known to lead to infrared divergences. This is the Linde argument [15,16] that diagrams with arbitrarily 
high numbers of loops all contribute to 0{g^). However it is possible to organize the 6th order contribution into the 
two terms cig^ and C2^^ In^ as it is done in lattice calculations [17]; then other resummation techniques [12] may give 
access to C2, while ci must be computed numerically. For a more extensive review of the resummation techniques in 
order to calculate quantities as e.g. Zqqb, we refer to [18]. 



D. From Two Loop Vacuum Diagrams to Free Energy at 0{g^) 

The gauge invariance of diagrams (3.18.B-E) for the vacuum case has already been discussed in Sec. II. Since the 
decomposition rule (2.6) has changed to (3.13), additional diagrams appear, namely 




and 




(3.20) 



Both diagrams are of 0{g^) and, as we will see in section V, cancel diagrams stemming from 0{g^) counter terms. 



Clearly these diagrams should be neglected in a calculation which is only accurate to 0{g^) 
result for Zqcb up to 0{g^) is a-independent. 



Hence the obtained 
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IV. ABOUT THE PROBLEM OF SHIFTING MOMENTA 



In the previous Sections we have shown how powerful the diagrammatic method of Cheng and Tsai is. However, 
we have not discussed one crucial technicality, which is the subject of the present Section. Namely if various diagrams 
are related to each other (like in (2.27)), shifts in the momenta have to be performed before the diagrams are really 
algebraically identical. This requires that the used propagators are not too singular. Otherwise shifting momenta 
might change the result of a loop integral and thus diagrams like (2. 28. A) and (2.28.B) would not cancel each other. 
Therefore naively using the Cheng-Tsai method and ignoring this problem can lead to serious inconsistencies. Indeed, 
one of the major results of the original Cheng/Tsai work was to point out the subtlety that the widely used naive 
versions of the different sets of temporal and Coulomb gauge Feynman rules are inconsistent with basic principles. 
This can be traced back to inappropriate shifts of loop momenta [4,9]. 

First we recall that momentum shifts are necessary to connect the gauge dependent parts of different diagrams 
(c.f. (2.27), (2.28)). As long as the involved propagators have no singularities shifting momenta in loop integrals 
causes no trouble. This might change if the propagators have poles. To see this we study the temporal gauge with 
its gauge condition 



^0 = 0. 

Naively one expects the ghost fields to decouple and the gluon propagator to be 

i ( .. W\ i k'y 



k"^ + IS 
D^\k) = D^^{k) =0 



^0 k'^ 



This is of course a special case of the general form (2.1) with the choice 



^0 



(4.1) 

(4.2) 
(4.3) 

(4.4) 



However the spatially longitudinal part of the free temporal propagator has a double pole ai ko = making 
momentum shifts difficult to control. To perform loop calculations one must introduce a pole prescription which 
obviously changes a^. We choose 



+ 



kn , 



ko + if] ko—ir]\ ^4 
Thus the longitudinal part of the gauge propagator becomes 



+ 



_{ko + if])'^ {ko - ir]Y _ 



(4.5) 



D,{k) :=^D,,{k) 



{ko + if])^ 
1 



^ 2ikoT] — rf 
k'^ + ie 



{ko - iriY 



1 + 



-2ikoT] — T] 



2 \ 1 



+ 



{ko + if])'^ {ko - if]) 



+ o{v) . 



(4.6) 



Note that this yields the principal value prescription of the /i^o-pole, which was frequently used for loop calculations 
in temporal gauge until it turned out that the results disagree with Feynman gauge calculation [19]. It reintroduces 
the ghost field 



G 



abc 



+i9f 



(4.7) 



and the temporal mode 



Doo = 



P +ie [ 2{ko+ir]y 2{ko - ir]^ _ 
Note that the expressions (4.6), (4.7) and (4.8) are easily calculated by inserting (4.5) in (2.1) and (2.8). 



(4., 
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Naively one expects that all the contributions involving ghosts and/or temporal gauge modes vanish for 77 — >- 
due to the appearance of rf in the numerators of (4.7) and (4.8). But this is not true as we will show now. Let us 
concentrate on a two-loop vacuum diagram with two longitudinal (L) and one temporal (O) propagator, as shown in 
fig. 4.9 [9]. 



k 




k + q 



If 77 is taken to zero before the loop integrations are performed there would be no contribution from this diagram 
since the temporal part (4.8) of the gluon propagator vanishes in this limit. But 77 serves as a regulator for the ko-po\e 
in (4.6) and should be taken to zero only after the energy integrations are performed. To check whether there is 
a contribution in this limit we restrict ourselves to the leading term in /y, i.e. we take into account only the first 
contribution of (4.6). Thus we have to calculate (c.f. 4.9) 

[ dko dpoik- p)og"' {k - p)o9"' ^DL{k) (p)Doo {k + p) 
J k^ 



^ J dko dpo {ko - 
1 

{k + py + ie 

By scaling ko = hioV^ Po = ttoT] we get 



Po)' 



1 



{ko + ir]y 
^2 



+ 



+ 



1 



{ko - ir]Y 
,2 



+ 



.(Po ^ir]Y {po -ir]Y _ 



{ko+Po+ir]Y {ko+Po-ir]Y_ 



dKo diTo {k.o - tto)^ 



+ 



1 



_{hio+iy {hio-iy_ 
1 



+ 



[{TTo+ir {TTo-iri 



+ 



^7y2(K: + 7r)2 _ (^ + p)2 [(^^ + 7^^ + ^)2 {^o + TTo 

1 1 " " 



— — ^ — / dhio diTo 

{k+p)^J 



{hio - TTo)^ 

1 



+ 



{hio + ^)^ (^0 - i) 



+ 



(tto + i)^ (tto - i) 



+ 



{ko + tto + i)^ {1^0 + TTo - i)^ 
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9 {k+p}^ 



(4.10) 



(4.11) 



which might be calculated with the contour method and does not vanish in contrast to the naive expectation that 
temporal modes should not contribute in calculations carried out in the temporal gauge. Whether temporal gauge 
modes and/or ghost fields yield non- vanishing contributions to observable quantities strongly depends on the pole 
prescription one uses for the ko-po\e in the longitudinal part of the gauge propagator. 

This shows that infrared singularities in the propagator must be handled with extreme care. In principle the same 
holds true for UV-singularities caused by propagators which do not vanish at infinity. A prominent example is the 
temporal mode of the Coulomb propagator [9] 



i^00(A:) = ^ ^ finite 
k 



for ko 00. 



(4.12) 



In the case at hand we should worry about infrared problems since they appear in the naive perturbation theory 
at finite temperature and cause the necessity of resummation. Indeed in our decomposition (c.f. (3.13)) 



™^ r<abdn \ 

k^m'^ + k^) " ^ 'k^ 



(4.13) 
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the terms at the r.h.s. have IR singularities. As we have shown in the previous sections the three terms are treated 
differently in further diagrammatic manipulations, e.g. the necessary momentum shifts for the first term of the r.h.s. 
of (4.13) is different from the ones for the third term. Thus loop integrals over the sum on the r.h.s. of (4.13) which 
are well-defined since the sum is IR safe have to be split in sums over loop integrals which might cause trouble. 

Fortunately these considerations are purely academic since the transverse projector ^^^^ — -^^1^^ in (4.13) is always 

accompanied by a factor k'^ which makes all terms at the r.h.s. of (4.13) IR save (c.f. (3.13)). Thus in our calculations 
performed in the previous sections no IR singular terms ever appear. Of course this holds only true for covariant 
a-gauges. Once axial gauges or the Coulomb gauge are used the considerations of this section might become relevant 
again. 

Now we turn to the three loop diagrams and prove the gauge invariance (more precisely the a-independence) of the 
QCD partition function as calculated up to 0{g^) using the resummation scheme outlined in subsection III A. 

V. GAUGE INVARIANCE OF Zqcd AT Oig"^) 

In this Section we prove the gauge invariance, more precisely the a-independence of the 0{g^) contributions to the 
partition function. Decomposing the corresponding diagrams of (3.18) with the help of the extended Cheng-Tsai rule 
(3.14), we obtain two different sets of diagrams: 

I. diagrams with the same graphical representation as the vacuum contributions but which have an implicit mass 
dependence stemming solely from the massive static gluon propagator D'p in (3.10). 

II. diagrams which contain explicit mass terms, i.e., they contain "crosses" or "cross-stars". These again split into 
two sets: 

(a) diagrams in which the crosses arise from the mass counterterms (3.12/3.19) to the three-loop diagrams. 

(b) diagrams in which the crosses arise from the additional terms in the identity (3.18) compared to (2.10), 
where an arrow pointing to a three-gluon vertex is decomposed. 

To check the gauge-invariance of the set of diagrams I, one follows closely that of the vacuum case. It is a 
straightforward application of the original Cheng-Tsai rules discussed in Sec. II. The most important technical steps 
are given in Appendix A. This leaves us with the gauge-invariance check of the set II. Here, we analyse in a first 
Subsection the structure of Ila and we find that the a-dependence of this set is given by the diagrams (5.2.E-G) below. 
In a second subsection, we turn then to the set lib showing that its 0{g^) contribution cancels exactly the remaining 
terms of Ila. 



A. Mass Counterterms to Three-Loop 

We discuss now the two-loop diagrams of the perturbative expansion (3.18) of Zqcd which occur as counterterms 
to the three-loop contribution: 
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Diagrams (5.1.C,E,F) contain the a-independent, static gluon propagator only. The a-independent part of the 
remaining two diagrams can be analysed with the help of (2.14) following the procedure of Sec. II C. Decomposing 
the a-dependent parts of diagram (5.1.B), we obtain: 




(B) (C) (D) 




This allows for the following cancellations: 

(i) Diagram (5.2.B) cancels the a-dependent part of (5.1.D). This is a consequence of the identity (2.17), relating 
a four-gluon vertex to two contracted three gluon vertices. 

(ii) The diagrams (5.2.C) and (5.2.D) cancel the a-dependent 0{g^) contributions (3.20) of the two-loop part 
(3.18.C-F) of Zqcd- That both sets of diagrams have exactly the same structure follows from the identity (3.15) 
which relates a contraction arrow on a cross to a cross-star. It remains to show that the prefactors of (5.2.C), 
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(5.2.D) and (3.20) allow for cancellation. 

To determine these prefactors, we first specifiy the symmetry factors attached to each diagram. Then, we 
symmetrize all diagrams according to the arguments of Sec. II D. For the three-loop contributions discussed 
here, this amounts to six gluon propagators to be symmetrized. Let us start with diagram (3.20). Their 
combinatorical factor is j2 and they represent an a-dependent part of the basic diagram (2.21.B). Symmetrizing 
the latter with respect to its six different gluon propagators gives another factor ^. The total prefactor of 
(3.20) is 4 • ^ • ^ — since always four differently symmetrized versions of (2.21.B) lead to the same version 
of (3.20) upon decomposition, e.g.. 




On the other hand, diagram (5.1.B) whose decomposition leads to (5.2.C) and (5.2.D), carries a combinatorical 
factor J. Symmetrization with respect to the possible gluon propagators gives a factor Consequently, the 
terms (5.2.C) and (5.2.D) have a total prefactor ^, too, and cancel the contribution from (3.20). We are left 
with the diagrams (5.2.E-G) which must be treated together with the set of diagrams lib. 

B. Three-Loop Diagrams with Explicit Mass Dependence 

The complete set lib of a-dependent three-loop contributions which have an explicit mass dependence is given in 
the following eqs. (5.4) and (5.5): 



G G G 




(I) (J) (K) (L) 
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(M) 



(N) 



(O) 



(P) 




(Q) (R) (S) 

All the diagrams in (5.4) are at least of 0{g^) as can be seen by splitting each loop integral into hard 0{T) and soft 
0{m) momenta. 

In contrast, some contributions from the "self-energy" -diagram (3.18.N) are 0{g^). In the following equation, all 
diagrams containing such 0(^^)-parts are labelled with letters. The remaining diagrams, being 0{g^) and higher are 
left unlabelled 
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diagram. For these combinations one can cancel all -parts with the help of two relations: 

(i) the definition of the Debye-mass (3.8/3.9) and 

(ii) the transversality of the one-loop self-energy. 

Let us first focus on (i). As an example we show the 0(^^)-parts of the diagrams (5.5.B/H/N/V) (note that the 
propagators with an explicit mass-dependence (the "cross-stars") are always soft): 




(A) (B) (C) (D) 

The definition of the one-loop Debye-mass (3.8/3.9) assures that the sum of the four diagrams above is just minus 
the diagram (5.2.E): 




Consequently, these five diagrams cancel. In the same way we can identify the 0(^^) -parts of the diagrams 
(5.5.A/G/M/W) with (5.2.F) and (5.5.C/I/0/X) with (5.2.G). Note that it is not surprising that the actual definition 
of the Debye mass enters at this stage of the calculation. If the mass m in the counterterm (3.12) is not defined properly 
via the static limit of the self energy, the result for the physical quantity of interest might be incorrect and gauge 
dependent. 

Now we turn to argument (ii). We group the diagrams which are still left together to four-diagram blocks so that 
the sum can be depicted as one diagram which has a one-loop self-energy insertion. This insertion is always Lorentz- 
contracted with A^, which contains k/^ for covariant gauges. Due to the known relation for the one-loop self-energy, 

A:^n^" = 0, (5.8) 
these blocks of four diagrams are zero. To give an example, we look at the diagrams (5.5.D/J/P/S): 




The same holds true for (5.5.E/K/Q/T) and (5.5.F/L/R/S). In this way 0(^^)-contributions of gauge dependent 
parts of the three loop diagrams are cancelled. Thus we have proven that the partition function up to 0{g^) is 
a-independent for an arbitrary covariant gauge. Of course it would be interesting to extend the proof to 0{g^). Here, 
however, all diagrams of (5.4) and (5.5) contribute and so far we haven't found a way to prove their cancellation 
diagrammatically. 
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VI. SUMMARY AND CONCLUSIONS 



We briefly summarize the necessary steps to prove the a-independence of a sum of loop diagrams in a general 
covariant gauge: 

• Start with the diagram which has only three-gluon vertices (like (3.18.C) for the two-loop and (3.18.F) for the 
three-loop case). Decompose one gluon propagator according to (3.10). Then treat the rest of the diagrams in 
the same way. 

• Use (3.14) and the equivalent of (2.12) as long as there are /i:^-arrows pointing at a three-gluon vertex, (2.16) 
as long as there are such arrows pointing at a fermion-gluon vertex. 

• Use (2.14) to pull arrows out of closed ghost loops (c.f. the last equality in (2.25)). 

• At this stage all a^-dependent diagrams which have neither contraction arrows nor four-point vertices should 
drop out. 

• Next, diagrams with different topology are related by the rules of Sec. II. To get an idea which diagrams are 
related it is instructive to contract the lines with contraction arrows. 

• In all remaining diagrams the momentum integrations have to be split up in hard and soft modes to flgure out 
which diagrams contribute to the order of g one is interested in. The current deflnition of the Debye mass and 
the transversality of the one-loop self-energy has to be taken into account. 

• All diagrams which have no counterparts so far ought to be of higher order in the coupling constant. 

• Finally symmetry factors have to be counted as outlined in Sees. II D and V A. 

In this article we have presented a method to check whether a resummation scheme preserves gauge invariance. This 
method is diagrammatic and relates systematically diagrams with different topology. This avoids tedious analytical 
calculations since it can be seen already at the diagrammatic level which algebraic expressions cancel each other. (This 
means that in a calculation the integrands cancels before performing the integrations.) Of course this diagrammatic 
method is less elegant than formal proofs of gauge invariance using Ward identities. On the other hand formal 
proofs are sometimes not straightforward or ambiguous due to singularities. Originally this diagrammatic method 
was invented by Cheng and Tsai [4] to check the consistency of different sets of Feynman rules for vacuum QCD. 
Here we have given an example that this method can be extended to check gauge invariance of physical quantities 
calculated in improved perturbation schemes. We have applied the diagrammatic method to the calculation of the 
free energy of flnite temperature QCD up to 0{g^). 

Strictly speaking we have only shown that the result for the partition function/free energy is independent of a 
for arbitrary covariant a-gauges. In principle it should also be checked that the result remains the same if one 
performs the calculation in a different class of gauges like e.g. axial or Coulomb gauge. Thus from a rigorous point of 
view further investigations are required to prove the gauge invariance of ^qcd- Of course the diagrammatic method 
outlined here is not restricted to covariant gauges but can be used for arbitrary gauges [4,20]. Note that then one 
has to deal with the momentum shift problems discussed in Sec. IV. In practice however a-independence is often 
used synonymously for gauge invariance. Indeed e.g.. the famous "plasmon puzzle" [1], i.e., the fact that the plasmon 
damping rate calculated without resummation turns out to be gauge dependent, already shows up for different values 
of a. 

We expect that the method presented here can be extended to more complicated resummation schemes (like e.g. "hard 
thermal loops" [2] where also vertex resummation is taken into account) and to electro-weak gauge theory (where the 
partition function in the vicinity of the phase transition is of special interest [21]). For the latter case we also refer 
to [22] where a similar diagrammatic method is used to minimize the computational effort for calculating S-matrix 
elements in vacuum quantum fleld theory. 
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APPENDIX A: SOME DIAGRAMMATIC S FOR O(p^) - GAUGE INVARIANCE 



In this Appendix, we give further details of how to prove the gauge invariance of the " vacuum hke " three loop 
contributions to the QCD partition function. 

The explicit mass dependent diagrams are discussed in Sec.V. In addition, the straightforward application of the 
diagrammatic rules of Sec. II to (3.18) results in approximately 100 a^-dependent diagrams. We have checked that 
the sum of these contributions cancels. For lack of space, we restrict our presentation here to a special example: The 
cancellation of all a^-dependent contributions to star diagrams. 

Cancellation of the Star 

Instead of presenting all diagrams we demonstrate for the star diagrams how the Cheng-Tsai rules work for three- 
loop diagrams. Then we outline briefly how all the other diagrams drop out. 

1. The star diagrams allow the following decompositions: 
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(A) (B) (C) 



Note that to emphasize the structure of the decomposition (2.1), Cheng and Tsai have introduced shghtly 
different notations for the two A-dependent terms, one coming with A{—k) and an overah + sign, the other 
with A{k) and an overah — sign. To ease a purely diagrammatic analysis, we depart here from this notational 
convention, we use A{k) = —A{—k) such that all A in the diagrams (Al) to (A5) are understood to represent 
factors A{—k). As a consequence, if the propagator (2.1) connects two vertices of the same type, then the 
orientation of the arrow does not matter, the two A-dependent contributions are equal and can be summed 
up. This leads to the factors 2 in front of (Al.C), (A3.C) and (A5.C), while in (A2) and (A4) respectively, the 
arrows point on vertices of different types and both A-dependent terms have to be dealt with seperatly. 
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Some diagrams already cancel on this level of the calculation. To be specific: 

(1.6.J) cancels (A3.C), (1.6.H) cancels (1.7.E), (1.6.P) cancels (1.7.C), (1.6.Q) cancels (A2.C). 

The decomposition and cancellation of the star diagrams containing fermions will be discussed below. 

3. The remaining star diagrams do not cancel each other but are related to diagrams with different topology. To 
make their cancellation clear we use the identities from section II B, specially (2.17), (2.18) and (2.19). 
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We start with relation (2.17) and place the A at the upper line near the left vertex. We connect the upper two 
lines with a three-gluon- vertex, we do the same with the lower ones, and we connect the remaining free leg of 
each vertex. In this way we get: 




With this relation diagram (1.6.B) cancels a-dependent parts from the diagrams (3.18.M) and (3.18.N) respec- 
tively. Now closing (2.17) in different ways and using a ghost-gluon vertex as well as a three-gluon vertex, one 
can find 12 more relations. With the help of these relations the diagrams (1.6.C), (1.6.D), (1.6.E), (1.6.L) are 
also cancelled by other diagrams. 

For the star diagram, we need no identity derived from (2.18). However it is necessary to cancel diagrams 
originating from (3.18.M) and (3.18.N). 

Finally we use (2.19). In an analogous way as above we get e.g. the following relation: 




(A) (B) (C) 

At all one gets five relations out of (2.19) which cause the diagrams (1.6.F), (1.6.G), (1.6.1), (1.6.K), (1.6.M), 
(1.6.N), (1.6.0), (1.6.R), (1.7.B), (1.7.D) to cancel each other. 

4. Now we discuss the decomposition and cancellation of the star diagrams with fermions. 

We only give a brief example of the cancellations. Decomposing e.g. (A4.C) in the way described above we get: 




(1.10) 



(A) (B) (C) 

On the other hand we close the fermion identity (2.20) to find 




(1.11) 



(A) (B) (C) 
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Hence the diagram (1.10. A) is cancelled by a diagram originating from (3.18.N). Similarly all a^-dependent 
diagrams containing fermions cancel. 

In the same way all a^-dependent diagrams of (3.18) drop out of the calculation of the partition function. 

APPENDIX B: FEYNMAN RULES 

We use the following Feynman rules for covariant a-gauges in Euklidien space-time: 



fermion propagator 



gluon propagator 



ghost propagator 



fermion-gluon vertex 




/ia 

ghost-gluon vertex b • * * * • • c F^^^ = —igf^^^k^ 




four-gluon vertex 




~\ahcd 
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In static resummed theory discussed in this paper, the gluon propagator changes to 



9[iv — Si_ioSi,oSkoO 



(l-a) 



(Bl) 



and one gets an additional mass-counterterm vertex: 



mass-counterterm vertex 



From these the vacuum rules are easily recovered: 

1. Replace the Matsubara sum by an integral as described in Sec. Ill A. 

2. Multiply every vertex with a factor i and every propagator with (— i). 

3. Correct the ^--factor in the ghost and gluon propagator to ^2^, in the fermion-propagator analogously to 



^—m-\-ie 



APPENDIX C: ADDITIONAL FERMION VERTEX 

Here we define the four point gluon-fermion vertex used in (2.20). We consider the following diagram, using the 
fermion rule (2.16): 

7c i 




m 
n 




=a 



(Cl) 



a\jk 



— mf-\-i£ ;^ -\- ^ — rrif -\- is J _ 



u,kl 



',kl' 



(C2) 



From eq. (C2) we can read off the four point fermion-gluon vertex: 



(C3) 



Note that • • ^ Tf/ 



28 



On sabbatical leave from Institut fiir Theoretische Physik, Universitat Regensburg, D-93040 Regensburg, Germany. 

E. Braaten and R.D. Pisarski, Phys. Rev. D42 (1990) 2156. 

E. Braaten and R.D. Pisarski, Nucl. Phys. B337 (1990) 569; B339 (1990) 310. 

R. Kobes, G. Kunstatter, and A. Rebhan, Phys. Rev. Lett. 64 (1990) 2992; Nucl. Phys. B355 (1991) 1. 

H. Cheng and E.-C. Tsai, MIT-preprint "The Theorem of Equivalence and The Anomalous Coulomb Interaction" (1986). 

Y.J. Feng and C.S. Lam, Phys. Rev. D53 (1996) 2115. 

P. Arnold and C. Zhai, Phys. Rev. D50 (1994) 7603. 

P. Arnold and C. Zhai, Phys. Rev. D51 (1995) 1906. 

C. Zhai and B. Kastening, Phys. Rev. D52 (1995) 7232. 

H. Cheng and E.-C. Tsai, Phys. Rev. Lett. 57 (1986) 511. 

H. Cheng and E.-C. Tsai, Chin. J. Phys. 25 (1987) 95. 

P. Arnold and O. Espinosa, Phys. Rev. D47 (1993) 3546. 

E. Braaten and A. Nieto, Phys. Rev. D51 (1995) 6990; Phys. Rev. D53 (1996) 3421. 
N.P. Landsman and Ch.G. van Weert, Phys. Rep. C145 (1987) 141. 

J. Kapusta, Nucl. Phys. B148 (1979) 461. 

A. Linde, Rep. Prog. Phys. 42 (1979) 389; Phys. Lett. B96 (1980) 289. 

J. Kapusta, Finite- Temperature Field Theory, Cambridge University Press, Cambridge, England (1989). 

F. Karsch, M. Liitgemeier, A. Patkos, and J. Rank, hep-lat/9605031. 
M. Achhammer, M. Thoma, and U.A. Wiedemann, in preparation. 

S. Caracciolo, G. Curci, and P. Menotti, Phys. Lett. B113 (1982) 311. 

S. Leupold, TPR-96-19, Universitat Regensburg preprint, submitted to Phys. Rev. D 

W. Buchmiiller, Z. Fodor, and A. Hebecker, Nucl. Phys. B447 (1995) 317. 

Y.J. Feng and C.S. Lam, hep-ph/9608219. 



29 



